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Cutoff effects in meson spectral functions 
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We study the lattice spacing dependence of meson spectral functions calculated in quenched QCD with domain 
wall fermions as well as clover Wilson fermions in quenched and partially-quenched QCD. We conclude that for 
lattice spacing a < 3 GeV all excited states appearing in the spectral functions are lattice artifacts. 



1. Introduction 

Temporal meson correlation functions in Eu- 
clidean time can be related to meson spectral 
functions by analytic continuation. Furthermore 
using the Maximum Entropy Method one can ex- 
tract the spectral functions from meson correla- 
tors calculated on a finite lattice 0. The method 
was applied at zero [1*21 '^^^ non-zero tem- 
perature [5l4i5i6) . In such an analysis one usu- 
ally uses all time slices in the correlation function 
not just the long distance part. Therefore one is 
sensitive to short distance physics which may be 
affected by lattice artifacts. The lattice distor- 
tion of the correlators at short distances trans- 
lates into lattice artifacts in the spectral func- 
tions for large energies. So far only Wilson gauge 
action and Wilson (clover) ferniion action have 
been used to study meson spectral functions. The 
lattice artifacts seen in numerical simulations for 
this combination of actions turns out to be quite 
different from what is expected in the free theory 
0. The purpose of this paper is to study cutoff 
effects of meson spectral functions using different 
gauge and fermion actions. 
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Figure 1. The pseudo-scalar (top) and vector 
(bottom) spectral functions with DWF at a^^ = 
2.0 GeV. The vertical line and band show the 
mass of the ground and first excited states. 



2. Numerical results 

We consider the temporal correlators of meson 
currents 



G(r,r) = ^(J(f,r)J(0,0)) 



(1) 



with J = q{x,T)Tq{x,T) and F = 75,7^ for 
pseudo-scalar and vector channels. For temper- 
ature T the spectral function is related to the 
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meson correlator by the integral relation 

nf f _cosMr-l/(2T)) 

G(r,T)^y^ d.a(c.,T)^-j-^^^^.(2) 

Although this relation is valid a priori only in the 
continuum, it has been shown recently that in the 
limit of a non-interacting theory it holds on the 
lattice HI as well. 

We use quenched domain wall fermon corre- 
lators on DBW2 gauge configurations from the 
RBC collaboration. Furthermore, we calculate 
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tree-level clover fermion correlators on quenched 
1-loop tadpole-improved Symanzik lattices and 
on partially-quenched 2-t-l flavor Asqtad fermion 
lattices from the MILC collaboration. 

To extract the spectral function we use the 
Maximum Entropy Method. We also parametrize 
the spectral functions as a sum of delta functions 
a{ui) = FiS{uj'^ — mf) and determine the pa- 
rameters Fi and using constrained curve fit- 
ting techniques 'S' which allows for better system- 
atic control of the results. 

We start the discussion of our numerical results 
with the case of domain wall fermions. Here we 
performed calculations at three different lattice 
spacings, — 1.3, 2.0, 3.0 GeV and valance 
quark mass around the strange quark mass. For 
the coarsest lattice we see a two peak structure 
in the spectral functions with the second peak lo- 
cated around 2GeV. The positions of the peaks 
are reasonably reproduced using a two exponen- 
tial fit. In Fig. ^we show our results for finer 
lattices, a"^ = 2.0 GeV for the pseudo-scalar and 
vector spectral functions. We see a three peak 
structure, with the first excited state being close 
to the position of the physical excited states of the 
77 and (f> mesons. The third peak is likely to be a 
lattice artifact. We also performed a simple two 
exponential fit of the corresponding meson corre- 
lators. The results of the fit are shown as the ver- 
tical line and as a band. The ground state mass 
from the fit (statisitcal error ^ 1% not shown) is 
in reasonable agreement with the MEM analysis 
while for first excited state there is larger statis- 
tical uncertainty indicated by the band. To see 
whether or not the second and third peak corre- 
spond to physical resonances or lattice artifacts, 
we calculated the spectral functions at = 3.0 
GeV. The results are shown in Fig. [5]which shows 
that the ground state peak did not shift (apart 
from a small shift due to slightly larger quark 
mass) while the second and the third peaks have 
moved roughly as a~^. This indicates that these 
peaks are lattice artifacts. Similar results were 
found for the Wilson action |1I2| . In Ref. they 
were interpreted as bound states of Wilson dou- 
blers. For DWF, however, such an interpretation 
is more difficuilt to apply. For this lattice spac- 
ing we have also performed a constrained fit with 



three exponentials. The ground state mass shown 
as the vertical line in Fig. [3 was obtained from 
a simple uncostrained fit which was then used as 
a prior for the constrained fit. The result for the 
second and third peak from the constrained fit is 
shown as bands in Fig. [5] Reasonable agreement 
between the MEM and constrained curve fits is 
found. 
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Figure 2. The pseudo-scalar (top) and vector 
(bottom) spectral functions with DWF at = 
3.0 GeV. The vertical bands indicate the error 
on the peak location as determined from a con- 
strained fit 



One may wonder whether the second and third 
peaks are artifacts of quenching. Absent dynam- 
ical quarks, all hadrons are absolutely stable (i.e. 
have zero width) no matter how large their mass 
is. If dynamical quarks are present, hadrons may 
decay, and the higher the mass the larger the 
width. Therefore higher lying resonances merge 
into a continuum. To check this we performed 
calculations in quenched and partially-quenched 
QCD using the clover action. More precisely 
we used gauge configurations with Symanzik ac- 
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tion generated by the MILC collaboration as well 
the dynamical 2+1 Asqtad configurations. In 
both cases the lattice spacing was a « 0.09/m. 
We use the tree-level value for the clover coef- 
ficient, Csw = 1, and K — 0.14 which corre- 
sponds to mps/mv — 0.68. Results for the spec- 
tral functions are shown in Fig. As one can 
see, the spectral functions for quenched and par- 
tially quenched QCD are very similar, indicating 
that for quark masses around the strange quark 
mass effects due to quenching are quite small. 
In addition wc have also calculated meson spec- 
tral functions wit HYP smeared clover action in 
2-1-1 partially-quenched QCD, i.e. we have ap- 
plied HYP smearing on the dynamical configua- 
tions before calculating the meson correlators. 
We have performed calculations at several quark 
masses. From extrapolation of the pseudo-scalar 
mass squared to zero we determine the critical 
hopping parameter to be Kc = 0.12718(2) (statis- 
tical error only) . The purpose of calculation with 
HYP smearing was to see whether the smoothing 
of the gauge fields removes the artificial peaks, 
and the corresponding results are shown in Fig. 
O As one can see from the figure, the lattice arti- 
facts are reduced but not removed completely by 
HYP smearing. However, as a by-product we ob- 
serve that number of exceptional configurations 
as well as the number of CG in the fermions ma- 
trix inversions are greatly reduced. In fact we did 
not observe any exceptional configurations in a 
sample of 59 configurations down to quark masses 
~ ms/8 (nis is the strange quark mass). 
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Figure 3. Pseudo-scalar (top) and vector (bot- 
tom) spectral functions for clover action in 
quenched and partially-quenched QCD. 
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